Defect Turbulence in Inclined Layer Convection 
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We report experimental results on the defect turbulent state of undulation chaos in inclined layer 
convection of a fluid with Prandtl number (a 1. By measuring defect density and undulation 
wavenumber, we find that the onset of undulation chaos coincides with the theoretically predicted 
onset for stable, stationary undulations. At stronger driving, we observe a competition between 
ordered undulations and undulation chaos, suggesting bistability between a fixed-point attractor 
and spatiotemporal chaos. In the defect turbulent regime, we measured the defect creation, anni- 
hilation, entering, leaving, and rates. We show that entering and leaving rates through boundaries 
must be considered in order to describe the observed statistics. We derive a universal probability 
distribution function which agrees with the experimental findings. 

(a)e = 0.08 (b) e = 0.17 (c) e = 0.17 
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Weakly driven, dissipative pattern-forming systems of- 
ten exhibit the spatiotemporally chaotic state of defect 
turbulence, where the dynamics of a pattern is dominated 
by the perpetual nucleation, motion, and annihilation of 
point defects (or dislocations) |jj . Examples can be found 
within striped patterns in wind driven sand, electrocon- |||| 

vection in liquid crystals [|, nonlinear optics §1, fluid FIG - !• False color images of undulation chaos § at 

convection |,|, autocatalytic chemical reactions [|| , and e = 30 °- Red/blue are rising/falling fluid and green is 
Langmuir circulation in the oceans @. The hope is that low-amplitude convection. Uphill direction is at the top of 

,ij ■ r,i i-rr i. i. x, i, the page and the region shown is the subregion of the cell 

the dynamics of these very diftercnt systems can be char- , . , , . „ . . , r . , .,, 

, , . . . , i . i . i ii used in the analysis. Positive detects point uphill, negative 

actcrizcd by a universal description which is based only , . , , .... , ,, , .„, , , rri .,, 

J detects downhill and the pattern drifts downward Jl 01 with a 

on the defect dynamics. period ^ 150r„. 

A first description of defect turbulence was given by 
Gil et al. M for a spatiotemporally chaotic state of the 

complex Ginzburg-Landau equation (CGLE). They pos- tracking all defects in a finite area of the convection cell 
tulated that the nucleation rate for defect pairs is inde- we measured, for the first time, defect creation, annihi- 
pendent of the number of pairs M , and based on the lation, leaving, and entering rates for a defect turbulent 
topological nature of defects the annihilation rate is pro- state. The observed pair creation and annihilation rates 
portional to M 2 . Through detailed balance, they showed agree with the predictions g. In an experimental system 
that these assumptions lead to a squared Poisson distri- where periodic boundary conditions do not apply, pairs 
bution for the probability distribution function (PDF) of defects are no longer an appropriate description. To 
of M. They also found agreement with this PDF in describe the statistics correctly, single defect leaving and 
simulations of the CGLE with periodic boundary con- entering rates through the boundaries must be consid- 
ditions. Rehberg et al. [|| measured the PDF of M for ered. Our data show that the entering rates are approx- 
defect turbulence in electroconvection of nematic liquid imately independent of the number of positive/negative 
crystals and found it to be consistent with the predicted defects N± and the leaving rates are proportional to N± . 
squared Poisson distribution. Later, Ramazza et al. in- We derive a universal PDF which reduces to the earlier 
vestigated a defect turbulent state in optical patterns and predicted squared Poisson distribution || when bound- 
found that their data was not conclusive. To date, stud- ary effects are negligible. Our measurements agree with 
ies in both simulation and experiment have relied purely this new distribution. In addition, by measuring both 
on comparisons of the PDFs. The gain and loss rates, N± and the undulation wavenumber we determined that 
fundamental to the universal description of defect turbu- the onset of undulation chaos coincides with the theoret- 
lence, have not been measured. In addition, effects due ically predicted onset for stable, stationary undulations, 
to boundaries were not considered. At higher driving, we observe competition between re- 
in this Letter, we report experimental results on the gions of ordered undulations and undulation chaos, sug- 
defect turbulent state of undulation chaos in inclined gesting bistability between a fixed point attractor and 
layer convection of a fluid of Prandtl number s=s 1. By spatiotemporal chaos. 



When an inclined fluid layer is heated from below and 
cooled from above for small temperature differences AT 
the fluid layer experiences not only a linear temperature 
gradient (as in Rayleigh-Benard convection) but also a 
symmetry-breaking large scale shear flow, with fluid ris- 
ing at the warmer and falling at the cooler plate [O. 
Above a critical temperature difference AT C , a pattern 
of longitudinal rolls aligned with the uphill/downhill di- 
rection is observed for intermediate inclination angles. 
For e > 0.02, where e = AT/AT C - 1 , these longitudinal 
rolls are unstable to undulations [Ojll]. While it was 
predicted theoretically that a perfect pattern of undula- 
tions is stable [[13| , a defect turbulent state of undulation 
chaos was found experimentally (L4J. Recently, undula- 
tion chaos was also observed in numerical simulations of 
the Navier-Stokes equations when random initial condi- 
tions were used |13| , while with controlled initial condi- 
tions undulations were found to be stable. As shown in 
the snapshot in Fig. |fp,, undulation chaos is characterized 
by a pattern of disordered undulating rolls, punctuated 
by point defects which carry a topological charge of ±2tt 
Jig] . Defects are nucleated pairwise in regions of low con- 
vective amplitude (green), enter and leave the system at 
the boundaries, preferentially move at right angles to the 
rolls, and annihilate as pairs. For e > 0.10, the defect 
turbulent state appears intermittently with ordered re- 
gions of undulations, suggesting a competition between 
the fixed point attractor (undulations) and spatiotempo- 
ral chaos (undulation chaos). This is shown by the two 
snapshots in Fig. Rib, c. The dynamics of the state is best 
viewed in the movies available at Ref. M. 

The experiments were conducted in the same appara- 
tus as Ref. p4J for an inclination angle of 9 = 30° . The 
fluid was compressed CO2 at a pressure of (56.5 ± 0.01) 
bar regulated to ±0.005 bar with a mean temperature of 
(28 ± 0.05)°C regulated to ±0.0003°C. For these param- 
eters, convection appeared at AT C = (1.763 ± 0.005)°C. 
The planform of the convection pattern was observed by 
the usual shadowgraph technique Q. A cell of height 
d = (388 ± 2)/mi and dimensions 101<i x 50d was used, 
for which the vertical diffusion time was t v sa 1.3 sec. 
All data was measured in a homogeneous central subre- 
gion of size 32<i x 25<i |fl7| . To track defects, runs of 100 
shadowgraph images were taken at a rate of 3 frames/sec 
(w t„/4) with a 10 bit IK x IK digital CCD camera. For 
17 values of e between 0.04 and 0.22 this was repeated 
for 600 (500, 400) runs at e < 0.07 (0.08 < e < 0.10, 
e > 0.12). Runs at the same e were separated by at 
least IOOtv, for statistical independence. A total of 1.5 
terabytes of data were analyzed. Each value of e was 
reached via a quasistatic temperature increase in steps 
of 0.001°C, waiting 700t„ between steps for transients to 
die out. We also conducted a sequence of measurements 
with quasistatic temperature decreases to check for pos- 
sible hysteresis, which was not observed. 

To lowest order, the planform of undulation chaos can 
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FIG. 2. (a) Roll wavenumber q and (b) undulation 
wavenumber p (black) with relative amplitude ipi/tpo (g ra y), 
as a function of e. Values are taken at the peaks of the PDF, 
with bars representing the width of the distribution at one 
standard deviation. Data taken while increasing temperature 
are shown with symbol A and decreasing with V- Dashed 
line represents prediction for onset from JOl. 



be captured by the real part of the two dimensional field 
ty(x,y) with 

* = tp Q e i{qx+ll>) + #1 (e i{qx+py) + e l{qx ' py) ) +■■■ (1) 

where the wavenumbers q and p describe the stripe spac- 
ing and waviness, respectively. Both the wavenumbers 
and the amplitudes ipo an d ipi vary slowly in space and 
time. For a perfect undulation pattern <f> is constant ev- 
erywhere, while for undulation chaos 4*(x, y, t) approxi- 
mates the deformation of the pattern due to the defects. 
The ansatz (f>(x,y,t) — J2 c » arctan( ^~^' ) describes the 
observed pattern well for constant q,p, and V'i/V'Oi where 
(Xi, yi) are the locations of defects and Cj = ±1 are their 
corresponding topological charges. 

In Fig. @ the peak values for the distributions of p, q, 
and V'l/V'o are shown as a function of e E3]. The values 
for q change little with e, while the undulation wavenum- 
bers can be fit by p oc -y/e — e Ul shown as a solid black 
line. This fit provides a value for the onset of undulation 
chaos at e u = 0.017 ± 0.001, which is in good agreement 
with e u = 0.016 predicted theoretically for the onset of 
stable stationary undulations (r|] . 

In order to better characterize the generic properties 
of defect turbulence, we determined defect locations and 
the corresponding topological charges in each shadow- 
graph image. A time trace of the number of defects over 
2x lO 5 ^ at e = 0.13 (10 4 defects) revealed that the undu- 
lation chaos state is statistically stationary and neutral 
with respect to topological charge. Nonetheless, any sin- 
gle realization of the pattern does not necessarily satisfy 
N + = A_ due to the boundaries. From the list of defect 
locations in each frame, we determined the mean (N±) 



(N±) - (N±) 2 in the chosen sub- 



and the variance a± 
region, plotted in Fig. H as a function of e. Here and 
below, we report results only for 7V + ; equivalent results 
were obtained for N- . 

For e < 0.10, N±(e) is linearly increasing. Extrapola- 
tion to N + = gives €d — 0.025 ±0.002, which is close to 
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FIG. 3. Number of positive defects in the subregion as a 
function of e, with bars representing the width of the distri- 
bution at one standard deviation. 




FIG. 4. Probability of entering, and leaving, creation, and 
annihilation (per 0.33 sec) as a function of iV+ for several e. 
Lines are fits to Eqn. H. 



the onset of undulations (e u = 0.017 ± 0.001) determined 
above from wavenumber measurements. The difference 
may be attributable to the finite size of the system. For 
e > 0.10, N±(e) decreases in the bistability regime, while 
no such decrease is seen in p. 

Two predicted generic features of defect turbulence f|] 
are the constant creation rate and the annihilation rate 
quadratic in the number of defect pairs M. By connect- 
ing defect locations in adjacent frames into defect tracks, 
and then associating tracks of opposite sign to locate cre- 
ation and annihilation events, we were able to perform 
the first direct test of the theoretical postulates, as well 
as the entering and leaving rates to/from the observa- 
tion area. For a system where single defects can enter 
and leave through the boundaries a pair-based descrip- 
tion is no longer justified. Thus, we separately consider 
the PDFs for N±. 

The creation rate C was observed to be approximately 
independent of N± and the annihilation rate agreed with 
A oc N+N-, by extension from A oc M 2 for the topolog- 



ically neutral case. Fig. [| shows the four gain/loss rates, 
tabulated using positive defects and events only. To low- 
est order, the entering rate is independent of the number 
of defects already in the subregion of the cell, suggesting 
that both creation and entering are random processes. 
The leaving rate was found to increase in proportion to 
N± as expected, since unlike annihilations such events 
do not depend on the presence of two oppositely-charged 
defects. 

The observed rates are approximately given by 



E(N±) = E Q 
C(N ± ) = C 
L(N ± ) = L N ± 
A{N ± ) = A Nl 



(2) 



which consider defects of a single sign only. The fits 
corresponding to Eqn. |2| are shown in Fig. |l[ for three 
values of e. It is important to note that this is a sim- 
plification. C(N±) and E(N±) show a weak, negative 
slope: defects are slightly less likely to enter or be cre- 
ated if there are already defects present. We also ob- 
serve a weak dependence of these rates on the number of 
oppositely-charged defects: E(N + ) is slightly diminished 
and L(N+) is slightly elevated for N + > N-. Thus, to 
fully describe the system, we would need to consider rates 
and distributions for both N + and iV_ simultaneously. In 
order to make analytical progress we will assume below 
that Eqn. @ are a complete description. As will be shown 
below, these assumptions approximate the data well. 

Using Eqn. we can construct the PDF for N = N±. 
Assuming a stationary distribution, detailed balance re- 
quires \oss(N)V{N) = gain(A - 1)T(N - 1) to describe 
the probabilities at adjacent N [091. For the relevant 
rates, 



V(N) 



E(N) + C(N) 



V(N-l) 



(3) 



L(N) + A(N) 

where a = °^ and (3 = 4 s -. By performing the recur- 
sion and properly normalizing the distribution, wc find a 
modified Poisson distribution 



v*A N ) 



a 2 



+N 



l0(2y/a)T(l + + N)N\ 



(4) 



where Ip is the modified Bessel function. For (3 = 
(a = (N 2 )) this PDF reduces to the squared Poisson 
distribution of Ref. ||. In our experiments, /3 rs 3 |20|. 

Fig. ||b shows that the squared Poisson and Poisson 
distributions match the data poorly. Both the modified 
Poisson distribution and a distribution obtained by per- 
forming the recursion of Eqn. on the raw experimental 
data from Fig. agree well with the experimental data. 
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neutrality, a description based purely on a single topo- 
logical charge was sufficient to describe the observed be- 
havior (2fJ. 

We thank W. Pesch and O. Brausch for sharing the 
results of simulations and stability analysis which have 
enhanced our understanding of the bistability. We would 
like to thank the NSF for support under DMR-0072077. 



FIG. 5. (a) Mean number of positive defects vs. a\ for e 
between 0.04 and 0.10. Black symbols are experimental; gray 
symbols are from Eqns. H and H; dashed line is the analytical 
relation for squared Poisson distribution fl8|. (b) PDF for 
6 x 10 4 defects at e = 0.07, with four distributions shown 
for comparison. O are the experimental distribution, solid 
line is modified Poisson, dashed is squared Poisson, dotted is 
Poisson, and + are from the raw rates in Fig. H. 



The modified Poisson distribution also correctly captures 
the mean and width of the PDF. Using the already nor- 
malized modified Poisson distribution (Eqn. H), we cal- 
culate (TV) = E N V a A N ) an d (^ 2 ) = Y<N 2 V a ,p(N) 
and obtain the analytic expressions 



(N) = 
{N 2 ) = 



- WVa) 
; Ip(2y/a) 



(5) 
(G) 



Again, the above equations reduce to the corresponding 
equations for the squared Poisson distribution in the limit 
(3 = 0. In Fig. Ha the dependence of the mean on the 
variance is plotted for the raw data and the modified 
and squared Poisson distributions. The modified Poisson 
distribution (Eqn. |5|) shows excellent agreement, while 
the squared Poisson distribution clearly fails to describe 
the experimental data. 

In summary, we find that the state of undulation chaos, 
and not the fixed point attractor to ordered undulations, 
is selected above the predicted onset for undulations. At 
higher driving, we find an apparent competition between 
the two attractors. 

Our result for the iV-dependence of the defect creation, 
annihilation, leaving, and entering rates and the modified 
Poisson (Eqn. Q) should be generic to any defect turbu- 
lent system. We have found agreement with the theoreti- 
cal predictions for creation and annihilation rates lq] and 
have extended the analysis to the experimentally relevant 
case of a finite system without periodic boundary condi- 
tions where it is important to consider both positive and 
negative defects (rather than pairs). Although we ob- 
served weak fluctuations away from topological charge 
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